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Abstract 

We generalize a recent model-independent form factor parameterization de- 
rived from rigorous dispersion relations to include constraints from data in 
the timelike region. These constraints dictate the convergence properties of 
the parameterization and appear as sum rules on the parameters. We fur- 
ther develop a new parameterization that takes into account finiteness and 
asymptotic conditions on the form factor, and use it to fit to the elastic tt elec- 
tromagnetic form factor. We find that the existing world sample of timelike 
data gives only loose bounds on the form factor in the spacelike region, but 
explain how the acquisition of additional timelike data or fits to other form 
factors are expected to give much better results. The same parameterization 
is seen to fit spacelike data extremely well. 

ll.55.Fv, 13.40.Gp, 25.80.Dj 



Typeset using REVT[TiX 



* Electronic address: buck@jlab.org 
^Electronic address: lebed@jlab.org 



1 



I. INTRODUCTION 



Dispersion relations in field theory represent nothing more than the application of 
Cauchy's theorem to Green functions, and yet information obtained from these identities 
continues to provide fascinating new insights into systems difficult or impossible to probe 
with other techniques. The natural advantage of the dispersive approach within QCD is that 
its incorporation of nonperturbative information is fully rigorous and model independent, 
since it works directly with hadronic Green functions in the contour integral. Thanks to 
Cauchy's theorem, the same quantity may be evaluated at a kinematic point where pertur- 
bative QCD provides the best physical description of the Green function. 

And yet the physical input is minimal. Nothing has been included from the fundamen- 
tal theory except for the value of the Green function at a single point and the existence 
of hadronic bound states of given quantum numbers. Clearly, QCD dispersion relations 
provide only a bare framework of conditions that Green functions must satisfy in order to 
be consistent with QCD (or any well defined field theory), namely unitarity in the form of 
quark-hadron duality, crossing symmetry, and analyticity outside of poles and cuts with lo- 
cations dictated by physical particle thresholds. There still remains a large space of possible 
hadronic Green functions, each of which satisfies the dispersive constraints. Nevertheless, 
the constraints can be surprisingly restrictive, especially when a small amount of additional 
physical input, such as data or the result of some model calculation, is included. It follows 
that constructing models manifestly satisfying dispersive constraints is an economical way 
of enforcing consistency with the basic features of QCD. 

The particular formulation of dispersion relations presented here focuses upon obtaining 
bounds on both the size and shape of electroweak form factors; this is a very old game, 
dating back into the early 1960s [ij, long before the advent of QCD. The essential notion 
is the computation of the current two-point correlator in two ways: at a point deep in the 
Euclidean region of momentum transfer t, and as an integral over the cut generated by on- 
shell hadronic states produced by the current. Components of the correlator are chosen so 
that both sides are positive definite, implying that the neglect of contributions from some of 
the hadronic states leads to a strict inequality on the hadronic amplitudes, i.e., form factors. 

The inequality thus obtained is expressed as a weighted integral of squared form factors 
over the kinematic cut in t corresponding to production of these hadrons. Assuming only 
that the branch point t+ of this cut is the lowest real value of t where the Green function 
exhibits non-analytic behavior, analytic continuation to the rest of the complex t plane 
carries the inequality to other kinematic regions. The inclusion of perturbative QCD Green 
functions in these studies was first made in Ref. 0. Subsequently, it was observed that 
the most general space of analytic functions satisfying the dispersive bounds obeys a simple 
parameterization 0, which led to a number of applications in the study of semileptonic 
decays, in which it is the weak current that appears in the correlator. 

Recent applications use only data below to restrict possible form factors. Surely 
data from t > must provide additional restrictions; however, for the semileptonic case 
this would require data from pair production via a weak current, which is beyond current 
experimental capabilities. Instead, we focus here on electromagnetic form factors, where 
data in the timelike region is abundant. Nevertheless, the mathematical techniques are the 
same. We show that data from the timelike region applied to the parameterization of Ref. 
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leads to a number of new sum rules on the parameters, which may be incorporated in an 
improved parameterization. These sum rules permit a more complete characterization of 
the form factor at all values of t. 

As an explicit example, we consider the elastic vr form factor and show that fits to the 
world sample of timelike data lead to predictions on the spacelike behavior of the form 
factor, particularly its normalization at t = and the pion charge radius. These predictions 
turn out to be very loose, but we shall see that this is the result of gaps in the data 
resulting from the fact little data has been taken near the pair production threshold t = 4m^. 
The same parameterization using spacelike data yields an excellent fit, in particular to the 
normalization of the form factor at t = and the pion charge radius. 

In Sec. II we review the derivation of the parameterization from dispersive bounds, with 
an eye toward its application to electromagnetic form factors. Section III demonstrates the 
extension of the dispersive bound to the inclusion of timelike data, and its effect on the 
parameterization. In Sec. IV we discuss issues of convergence of the parameterization in 
the timelike region, and observations made in this inquiry lead us, in Sec. V, to develop 
an improved parameterization. Section VI presents the result of fits to the elastic vr form 
factor and discusses at length the quality of the fits and extrapolations from timelike to 
spacelike regions. We also compare the results to direct fits to spacelike data using the 
same parameterization, as well as fits from the literature. Section VII summarizes possible 
improvements and concludes. 



II. REVIEW OF THE DISPERSIVE APPROACH 

Much of this discussion is patterned on that in Ref. |^ and references therein, but appears 
here again for clarity in understanding what is to follow. We begin with the QCD two-point 
function of a vectorlike current J'^; in the electromagnetic case, this is simply the conserved 
vector current qj^q. The polarization tensor is defined by 



n7(g) = ^iq'q"" - gV'^)nJ(g2) + ^n^(g^) = z / Ae^<^-(0|TJ'^(x)jt^(0)|0), (2.1) 



and in the case of a conserved current J, the polarization function 11^ vanishes. The re- 
maining function 11^ is rendered finite in QCD by making two subtractions, leading to the 
dispersion relation 



Xj(r) = oTJT^ = - / dt- (2.2) 

2(9(g2)2 7r Jo {t-q^Y ^ ' 

The function Xjil"^) be computed reliably in the perturbative QCD operator product 
expansion for values of far from the kinematic region where the current J can create 
on-shell hadronic states; for any light quark current and = — g^, this condition reads 

Inserting a complete set of states X into the two-point function relates IlJ to the pro- 
duction rate of hadrons from a virtual photon. 



ImnJ = lB2^)V(g-px)|(0|J|X)|% (2.3) 

X 
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where the sum is over all hadronic states X with the same quantum numbers as the current J, 
weighted by phase space. It follows from the dispersion relation (p.2|) that the perturbatively 
evaluated Xj{q^) equals the integrated production rate of 7* — X weighted with a smooth 
function of momentum transfer squared t. Since the sum is positive semidefinite, one may 
restrict attention to a subset of hadronic states to obtain a strict inequality. In the case 
discussed in this paper, J'^ = ^U'y'^u — ^d'y^d, and we restrict to X = 7r+7r~.[] This places an 
upper bound on the electromagnetic vr form factor F{t) in the pair-production region that 
takes the form 



1 



(t-q 



2A3 



(2.4) 



from the dispersion relation Eq. ( |2.2| ). Here W{t) is a computable function of t that depends 
on phase space and the quantum numbers of the particular form factor under consideration. 

Using analyticity to turn ( p.4[ ) into a constraint in the spacelike region in t requires 
that the integrand is analytic below the pair production threshold t < t+. To do this, we 
introduce a parameter tg < t+ and a function 



Z{t]ts 



(2.5) 



that is real for t < t+, zero at t = is-, and a pure phase for t > t+. Any poles in the 
integrand of Eq. ( |2.4| ) can be removed by multiplying the integrand by various powers of 
z{t]ts), provided the positions of the sub-threshold poles in F{t) are known. Each pole 
has a distinct value of t^, and the product z{t; t si) z{t; ts2) ■ ■ ■ serves to remove all of them. 
Such poles arise as the contribution of resonances with masses below to the form factor 
F{t), as well as singularities in the kinematic part of the integrand. After determining these 
positions from the hadronic mass spectrum, the upper bound on F{t) becomes 



POO 

/ dt 


dz(t; to) 


dt 


Jt+ 



m;to)Pit)Fit)\'<l 



(2.6) 



where the weight function 0(t; to) (known as an outer function in complex analysis) is given 
by 



0(t;to) = P(t) 



W{t) 



\dz{t;to)/dt\x^{q'){t~q'r 



(2.7) 



The factor P(t) is a product of z(t;ts)'s and yz(t;ts)'s, with tg chosen to remove the sub- 
threshold singularities and cuts in the kinematic part of the integrand, while the Blaschke 
factor P(t) is a product of z(t;tp)^s with tp chosen to be the positions of sub-threshold 
poles in F(t). The functions 0(t; to) and P(t) also depend on g^, which we leave implicit for 
notational simplicity, while to is the (yet to be chosen) value of t for which z(t; to) = 0. 
The dispersion inequality expressed in terms of z reads 



^The neutral pions do not appear here because of charge conjugation. 
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-i. / -\<Piz)Piz)Fiz)\'<l, (2.8) 
Zm Jc z 

where F{z) means F[z(t;to)] and so on. Lacking poles, the quantity (f)(t;to)P{t)F(t), is 
expected to be analytic within the unit discj^ and may be expanded in a set of orthonormal 
functions that are simply powers of z{t; to). The result of expanding in z(t; to) is an expression 
for F{t) valid even in the spacelike region, 

-I oo 

^(^)= P(,^^(, ,^ T.^nz{t;to)- , (2.9) 
P{m{t;to) „=o 

where, as a result of Eq. (|2.8|) , the coefficients a„ are unknown constants obeying 

|a„|2 < 1 . (2.10) 



oo 
n=0 



The functions P(t), (j)(t;to), and z{t;to) are real by construction for t < t+. Moreover, 
physical cuts in the form factor, which lie on the real axis, generate discontinuities only in 
the imaginary part of F. Since there is no physical distinction between the upper and lower 
complex half-planes, one must have \F[t+ie)\ = \F{t — ie) \ for all real t, and consequently the 
form factor satisfies the Schwarz refiection principle.^ It follows from analytic continuation 
away from the cut that the form factor is real on the real t axis below threshold, which maps 
to the real z axis, and therefore that the coefficients a„ are real. 



III. CONSTRAINTS FROM THE TIMELIKE REGION 



The constraints derived in the previous section require no input from the value of the 
form factor F{t) for t > t+. Clearly, if such information is supplied from another source such 
as direct measurement, the bound provided by, e.g., ( |2.4| ) should be strengthened. In terms 
of the equivalent bound of Eq. (|2.8| ), such data appears on a segment of the unit circle in z. 
Suppose that one is given data for F(t) from threshold up to some tu. According to (^.51), 
the segment (t+ =F ie, tu =F ^e) written in terms of angles from z = e*^ occupies the circular 
segments 9 G (6'„,7r) and (— vr, —9u), where 



cos 



-1 



{tu-t+) - {t+-to) 



t,, — tn 



(3.1) 



^Of course, subthreshold singularities due to multiparticle or anomalous thresholds must be con- 
sidered. In many cases [Q-^ these singularities treated as cuts are numerically unimportant, and 
in the Kit Q, tt+tt^, or KK cases they are absent. Moreover, we here correct an oversight of these 
previous works: For on-shell multiparticle resonances of the correlator or triangle- type anomalous 
threshold diagrams, the kinematics of the internal particles is completely fixed, meaning that the 
singularity in t is a pole, not a cut. 

^Alternately, the form factor is initially defined only in the upper t half plane, and the Schwarz 
reflection principle continues it into the lower half plane. 
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The additional information tells us about the integrand of Eq. ( p.8| ) directly on the unit 
circle; here the Blaschke factor P{z) is unimodular, and so we may express this input as 
follows: 



n=0 



\n=0 

oo 



Vr?i=0 



oo 



E «n E «m^' 
n=0 m=0 

oo oo 



E E '^ni [cos(n — m)6u + i sm.{n — m)6^ 



(3.2) 



n=0 



m=0 



Since the value of the modulus \(f){z)F{z)\ is invariant under reflection about the cut (which 
corresponds to du —>■ —Ou), the left hand side of the corresponding expression for z = e~*^" is 
identical, while the right hand side flips the sign of the sine term; this term must therefore 
vanish. Equivalently, noting that the left hand side is real and recalling from Sec. II that the 
ct^ s are real, the purely imaginary sine term must again vanish. Yet another way of seeing 
the same result is by using sin(n — m)9u = smn9u cos m9u — cosn9u smm9u and writing the 
full term as 



an smn9u E '^m cos m6'„ ] - [Y ^nCosn9u \ i ^ sin m6'„ 



(3.3) 



\n=0 



\m=0 



\n=0 



\m=0 



Z 



which is pure imaginary and vanishes trivially if all the real. We are therefore left 

with one nontrivial sum rule: 

2 oo oo 

^''-)Fiz = e''-)\ = 5]a„ ^ a„cos(n-m)^„, (3.4) 

n=0 m=0 

We arrive at the mathematical result that allows for anaylsis of timelike data. Its deriva- 
tion has been exceptionally straightforward, but a number of comments are in order before 
proceeding. 

First, the derivation is completely consistent with the constraint (|2.1CI|) . Indeed, inte- 
grating cos(m — n)9 over the unit circle gives 2'7i6mn, and thus ( p. 41) integrates to ( ^.10| ), 
using (p?8|). 

Second, although the first line of ( |3.2| ) is written as a perfect square, exploiting the reality 
of the Ori's means that it is more convenient to combine the two sums as in ( p.4|) . This is 
because the sum rule applies to the magnitude of F for t > t+. Although phase data on 
above-threshold form factors exists in some cases, typically it is data for \F\ above threshold 
that is presented in experimental papers. Indeed, the very nature of our dispersion relations 
[see ( p.4| )] precludes our use of such phase information. 

Third, we immediately see from taking particular values of 9u in ( p.2|) , or even from its 
first line, two very interesting special cases: 

oo 
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TT 



{t = U) : 
9u = {t^ +oo) : 



(+1)^1)1 



E«n(-ir 

n=0 



E«" 

n=0 



(3.5) 
(3.6) 
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Clearly, conditions on or F at these special points provide strong constraints on the 
parameters. We explore this next. 



IV. CONVERGENCE OF THE PARAMETERIZATION 

In the last section we manipulated infinite series without regard to the fine points of 
convergence. Does the derivation of our sum rules suffer from the possibility that some of 
these actions are ill-defined? In this section we show that this is not the case. 

The central convergence problem may be posed as follows: In deriving (|2.10|) , we use 
the analyticity of the combination P{z)(j){z)F{z) inside the unit circle to expand it in the 
Taylor series J2n ^n^"- Since physical information "along the cut" actually appears at values 
of momentum transfer an infinitesimal distance from the genuine cut [t ± ie] , it follows that 
this data is actually infinitesimally inside the unit circle in z, where the geometric series 
X]„ z"' is still (barely) convergent. The finiteness of analytic complex functions on a compact 
region such as the disc \z\ < 1 — e, e > 0, gives one sense in which the sum J2n(^nZ"' has 
meaning as a finite number, through geometric convergence. On the other hand, for |z| = 1 
the expression J2n(^nz'^ is defined as the value of the quantity P{z)(j){z)F{z) along (the 
appropriate side of) the cut, which exists since the only place where this combination is ill 
defined, the cut t >t+, has been mapped so that its two sides are sent to the two separated 
halves of the unit circle in z. These two expressions for J2n(^nZ^ at \z\ = 1 are equal due 
to a theorem of Abel's: In words, since analyticity demands that the power series X^n ^n^" 
converges for all \z\ < 1, and since its value at 1^1 = 1 is defined, the \z\ —>■ 1 limit of 
the former equals the latter. Therefore, the expansion as an infinite power series in z with 
1^1 = 1 is true as a formal statement, but does it have meaning as a useful series with rapidly 
converging partial sums when we do not invoke the technicality that we are "just inside" 
the unit circle? 

The utility of the expansion in z, when continued to the sub-threshold region for semilep- 
tonic decays, depends on its geometric convergence, using l-z] < 1 inside the unit circle and 
the boundedness of the a„'s via ( ^AOl) . It follows in that case that only the first few a„'s 



are required to describe the form factor over the entire semileptonic region However, 
\z\ = 1 on the unit circle, and so geometric convergence fails completely. A priori it seems 
that an arbitrarily large number of anS is required to describe the form factor in the region 
t > t+. 

Yet things are not so bleak. We now show that the analyticity structure plus some mild 
physical assumptions tells us much about the convergence properties of the parameterization, 
even on the unit circle \z\ = 1. 

To begin with, in order to use Parseval's theorem rigorously to prove (|2.10|) , one must 



be able to exchange the order of sum Y.n ^n-z" and integral in z; for this purpose, we de- 
mand the sufficient condition that X^n'^n-^" is uniformly convergent in z for \z\ = 1. Now, 
\P{z)(j){z)F{z)\ is by physical assumption a bounded, smooth, continuous function on the 
compact region represented by either half of the unit circle.^ However, even with these 



^The functions \(l){z)\ and |-P(2;)| = 1 have no singularities on \z\ = 1, as discussed below, so our 



7 



assumptions, uniform convergence could be spoiled by a region of z where the form factor 
is not smooth in z even though it is smooth in the original kinematic variable t. The only 
region fitting this description is t — +cx), which is compressed into the finite region z +1, 
\z\ = 1. That is, a form factor could continue to oscillate smoothly in t all the way out to 
infinity, but mapped to z this behavior would appear to oscillate wildly. We therefore make 
the additional physical assumption that becomes featurelessly smooth for t +00 

in order to guarantee uniform convergence in z and thus the proof of ( |2.1CI| ). 

Now we may consider the full power of ( p.lO| ). Rather than using only the boundedness 



of each term as in the geometric case, we recognize that the absolute convergence of the sum 
of |a„p implies that for sufficiently large n, |a„| falls off faster than Notice that, 

although this new criterion restricts the pattern of a„'s, the condition of convergence of 
J2n ci'nz"' is still strougcr. For, one could imagine, for example, that the true form factor has 
a„ = +l/\/C(3/2)n^/^, which satisfies J2n'^n — ^ t)ut clearly leads to divergence of Y.n'^nZ^ 
at z = +1. To proceed further, we need additional input. 

Up to this point, we have ignored the specific form of the function to) . In the notation 
of [^, one begins with 

imn^ > -^(t - t+)t(t - rlF(t)p0(t - (4.1) 

where a, b, c, and K are integers specific to the form factor under consideration, and nj 
is an isospin Clebsch-Gordan factor. In the electromagnetic case for pseudoscalars, only 
the form factor analogous to f+ in B ^ D appears, for which a = 3, b = 3, c = 2, and 
K = 48. The kinematic factors are simply an expression of two-body phase space, with 
t± = (M ± m)^ for M ^ m decays, which for the elastic case simplifies to t+ = 4M^, 
t_ = 0. The manipulations described in Sec. II lead to 



nr t+ ~ t \ ^ f I I \ r N ^ 

IT- 1 / /J J . I J I \ /J j_\ 4 



iU-t + ^u) i^^u - 1 + + ■ (4.2) 

A more convenient expression for \(j){t]to)\ with t > t+ reads 



mt; to)\ = t-'/\t - - to)-'^'Vt^o (t + QY'^\ (4.3) 

Note in particular that (J) ^ (t — t^Y as t — > t+, while as t ^ +00. The vanishing 

of the first limit indicates powers of the spatial momentum [|p|^ for the vector form factor] 



statements refer to the behavior of |F(2;)|. Resonances above threshold have finite widths and 
residues, so |-F(2:)| is also smooth, and we assume that |-F(-z)| is finite as t — > +00, i.e., z +1. 
Finally, the one-sided limit z — > — 1 along the circle is assumed to exist. 

^Of course, such statements here and below refer to the magnitudes of a„ in a statistical sense. 
Individual terms may deviate above or below the given n dependence, but the overall pattern of 
the terms is bounded by the stated behaviors. 
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of the pair created near threshold, while the vanishing of the second limit indicates the 
asymptotic unitarizing behavior of |p|^/v^. Making only the mild physical assumptions 
that the form factor F is not infinite at threshold and does not grow as t +00 [already 
assumed in proving (|2.1(]| )], ( p.5|) and (|3.6| ) give the sum rules 



£a„ = 0, £a„(-lf = 0, 



n=0 



n=0 



or equivalently, 



0, 



n=0 n=0 



0. 



(4.4) 



(4.5) 



Further improvements are possible if one considers the particular nature of the vanishing 
of 0(t; to) as if: — i> t+. Since the parameters a„ are defined with respect to the variable z, 
this is most obvious if one considers as a function of z rather than t, in which one sees 
that (j){z) oc (1 + z)*^""*"^)/^. For all form factors studied in that paper, a = 1 or 3, which 
originates from the suppression of pair production by Ipl'* at threshold. For a = 1, not only 
does 0(z) vanish at t = t+, but 0'(z = —1) is finite, while for a = 3, (t)\z = —1) vanishes 
and (f)"{z = — 1) is finite. To be explicit, 0(z) for the vector form factor written in z assumes 
the form 



0(z;to) 



[l + z) 



5/4 



'l-^(l + .) + (l-.) 



n-1/2 



X 



{l-z) + Jl--^{l + z) 



(4.6) 



Writing J^n^^nZ"" = P{z)(f){z)F{z), we obtain 

00 

Ean(-1)" = (P0F)(-1), 

n=0 

00 

Y: a„n(-l)" = (P0F)'(-1), 

n=0 

00 

^ a„n(n -!)(-!)"= (P0F)"(-1). 



(4.7) 



n=0 



Since one may confirm that P, P', and P" are finite on the circle, making the physical 
assumption that the nth z derivative of F{z) near z = —1 is no more singular than (l + z)"" 
gives 



5:a„(-l)" = 0, 

n=0 
00 

5]a,n(-l)" = 0(a = 3); 

n=0 

00 

J2ann{n-l)i-ir ^00 {a = 3). 



^00 {a = 1), 



(4.8) 



n=0 
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The second statement improves upon the constraint of the previous sum rule by implying 
|a„| — > faster than 1/n, rather than 1/y/n: Now we may state that ^n'^nZ^ is absolutely 
convergent. Thus one may reorder terms in the sums, as was implicit in the proof of Eq. ( p. 41 ). 

Moreover, in the case of the vector form factor (a = 3), the second and third expressions 
together tell us that in fact a„ — * at least as fast as This leads to a tremendous 

improvement in our knowledge of convergence of X^n o-nZ^ on the circle.^ Where before we 
knew only that convergence occurred, now we can quantify how fast. The relative error of 
the partial sum Yl,n=i ^/ri^ from its exact value of C(2) = 7r^/6 is bounded by 6/(7r^iV). We 
find that only 3 terms are required for a 20% relative error, 6 for 10%, and 60 for 1%. In 
principle, it should be very a simple matter to fit the coarse structure of the form factor 
over its entire kinematic range with very few parameters; even such a structure as the large 
resonant p peak may be accommodated by this quasi-Fourier analysis. However, it turns 
out that this expectation is unfulfilled, for reasons that we now discuss. 

The problem is that we do not know which of the most important in this 

convergence. That is, they are not necessarily the first few {oq, cti, ■ ■ ■ , ctAr}, as for the 
semileptonic decays or more generally for any spacelike factor, where moments about z = 
[t = to) dominate. However, the constraint (|2.10| ) tells us that, if this set is appreciable in 
magnitude, there is little room for higher a„'s; this leads to a rapid convergence in n, and 
consequently, one obtains a unified picture of the form factor spanning both spacelike and 
timelike regions. The alternate possibility is that one could have many a„'s much smaller 
than the saturation paradigm 1/n^, so that convergence is painfully slow in n; in this case, 
the fit to timelike data voraciously demands ever higher a„'s and is relatively insensitive to 
the lowest coefficients, so that the spacelike data is relatively unconstrained by the timelike 
data. A detailed fit to data is required to determine which scenario is realized. 

We pause momentarily to review the assumptions made to obtain appropriate conver- 
gence properties of the parameterization. Convergence of the series J2n c^n^" for t > t+ is 
obtained by requiring to be a bounded, smooth, continuous function, while in obtain- 

ing the sum rule ( |2.10| ), we additionally assume that becomes featurelessly smooth 

as t ^ +00. Absolute convergence of the series uses the particular forms of P and 0, and 
the assumption that F'{z) near z = —1 is no more singular than (1 + z)"^. Finally, the 
convergence of the series like 1/n^ depends on the form of (f){z) for the vector form factor 
and the assumption that F"{z) near z = — 1 is no more singular than (1 + 2;)"^. 



V. A NEW PARAMETERIZATION 



After the detailed discussion of the series X)r 



a„z' 



and its properties, it may seem incon- 



gruous to introduce a new parameterization for the form factor. Yet it is entirely appropriate 



to do so, since it is more natural to incorporate the sum rules (4^) and (O) directly into 
the parameterization than to impose them by hand. The sort of trouble one might encounter 
by using the parameterization in a„'s becomes clear with refiection upon the content of the 



^A similar treatment with powers of (1 — z) about t = +oo leads to the same convergence for the 
scalar form factor in processes with non-conserved currents. 
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sum rules (^7^). One may impose explicit constraints on a truncated set {ao,ai, . . . ,aN} 
to enforce the vanishing of the first two of these as well as the first expression in (|4.4| ), for 
example by fixing oq, ai, and 02 by means of the sum rules and the values of {03, 04, . . . , aAr}. 
But then these sum rules, which came from certain limits of the parameterization in z, might 
be accomplished in perverse ways. For example, the final expression in ( |4.8[ ) may be finite 
but exceptionally large, leaving essentially unbounded. 

The cure for such phenomena is straightforward. One simply factors out the appropriate 
behavior determined by (j){z) from the parameterization: 



00 

E 

n=0 



+ 



a/2 



n=0 



(5.1) 



defining the series J^n^nZ"'- Since the original series I]„a„2;" and also the prefactor (1 + 
z)~^(l — z)~^/^ are analytic in z inside the unit circle, the same holds for X^n^n^"- The 
requirements on the singularity structure of F{z) as z — ±1 discussed in the previous 
section are carried verbatim to J2n ^nz"' since 1 / (j){z) has no singularities on the unit disc 
except those removed by the prefactor, while the sum rules ( [4.4|) , ( [4.8|) are automatically 
satisfied by the inclusion of the 2;-dependent prefactor. The parameters 6„, like a„, are real 
since they in particular describe the form factor in the spacelike region, where F{z) and z 
are real. On the other hand, the b^s no longer satisfy any particular constraints except that 
J2nbnZ'^ is analytic.^] The extra prefactor is just that appearing explicitly in 0(2), so one 
obtains the effective weight function 



[z;to)=<P{z)il + z)-\l-z)-'/' 



•l-f{l + z) + {l-z) 



-1/2 




so that 



mz)\ n 



The sum rule that relates timelike data points to the parameterization, Eq. 
assumes the parallel form 



00 00 



z = e''^-)F(z = e'^-' 



(5.2) 

(5.3) 
thus 

(5.4) 



n=0 m=0 



Nevertheless, it is \(l){z)F{z) \ and not \(t){z)F{z)\ that appears in the dispersion integral, and 



so the original sum rule (|2.10|) in the new basis is not just X^n b\<l, but assumes the more 
complicated form 



''^An additional constraint occurs if a particular t 
assumed. 



+00 or t — > t+ behavior of the form factor is 
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J2^ri^bm5{m,n) <1, (5.5) 



n=0 m=0 

where 



64 2 



5(m,n) = — (V -22p + 45) n [(2« + 1)^ -p] , (5.6) 

i=0 



with p = 4(m — n)^. In the original sum rule (p.lO|) , where the analogue of S{m,n) is just 



Smn, one clearly requires each |a„| < 1, and every nonzero a„ serves to limit the size of 
succeeding terms. However, it is not quite so simple to determine at a glance whether the 
parameters 6„ obtained from fitting to data or a model satisfy the original dispersion relation 
( [^.81 ). Moreover, without the sum rules that produced a strongly convergent series in a^, the 
convergence of the series in 6„ follows no guaranteed pattern but still may give an adequate 
fit. We will see that this is in fact true for the vr form factor. 

An important aspect of the transformation between the parameters a„ and 6„ represented 
by (|5.1|) is that both series require only a small number of parameters to describe the form 
factor in the vicinity of z = 0; such behavior is necessary for a useful and minimal description 
of the spacelike form factor around t = to- Specifically, one finds 

ao = bo, 

3 

ai = bi + -bo, 

3 1 

^2 = ^2 + ~ 4^0, (5.7) 

and so forth. The essential conundrum in obtaining a parameterization valid in both space- 
like and timelike regions is that, in the spacelike region, or at least near t = to where data is 
abundant, it cannot be very different from the parameterizations in or 6„ discussed above. 
That is, it must have a rapidly converging Taylor series about z = 0. On the other hand, it 
must be able to recognize strongly localized structures in the timelike region \z\ = 1, such 
as the p resonance in the present case. We return to this point in the next section, once we 
have exhibited our empirical results. 

To summarize, the parameter basis a„ admits a number of sum rules from timelike 
constraints which give rise to a strongly convergent behavior. It also admits the very simple 



dispersive bound ( p.lO|) , but may exhibit pathological behavior in fits to timelike data. The 
parameter basis 6„ has all of the sum rules built in and is much more stable in timelike fits, 
but its dispersive bound (|5.5| ) is much more complicated, and in general may converge much 
more slowly. Both bases are useful in the spacelike region, while it requires a fit to data to 
determine whether either is useful in the timelike region. 



VI. THE ELASTIC vr FORM FACTOR 

Our intent is to study the possibility of using the considerable tt~^7t~ production data 
in the timelike region t > t+, from which the form factor |F(t)p is extracted, to fit to the 
parameters bn as determined from the analytic structure of the QCD dispersion relation 
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( p.2| ) and expressed by (|5.2| )- ([5l^ ). Of particular interest is the number of parameters 
necessary to give a good accounting of the data, and what these parameters tell us about 
the shape of the form factor in the spacelike region, where there exists considerable data 
from TT"*" elastic scattering. 

A. Inputs 

Let us first consider the calculation of the current two-point correlator, as represented by 
its subtracted form XviQ"^) iii (|2.2| ). As noted in Sec. II, this quantity possesses a well defined 
operator product expansion in inverse powers of Q^. The first few terms are particularly 
simple in the case of the light quarks u and d. Neglecting subleading mass corrections, 

(6.1) 

Note that the qq condensate has been included, despite the vanishing of the quark mass m, to 
indicate its place in the expansion. The electromagnetic current bilinear J'^ = ^u'y^u — ^d'y^d 
induces an extra factor of (+2/3)^ + (—1/3)^ = 5/9 in Xv through the quark charges. Our 
goal is to choose as small as possible in order to maximize the stringency of the dispersive 
bound through perturbative QCD input, but large enough that the perturbative expansion 
remains valid. 

Although numerical estimates for the first few condensates certainly exist in the lit- 
erature, we use them principally to establish that region in for which convergence of 
the expansion is satisfactory and to obtain a numerical uncertainty on the lowest-order re- 
sult. To be specific, we use {asG^/n) = 0.02-0.06 GeV^ and the expression for as{Q^)/Ti 
from the three-loop beta function with n/ = 3 and A^g ^ = 380 ± 60 MeV. Values for the 
as correction and the gluon condensate relative to leading order are presented in Table I 
for various values of Q^. We learn that corrections become quite large for < 2 GeV^. 
Although the 0{al) correction and the gluon condensate have opposite signs, large uncer- 
tainties on the latter prevent one from knowing how complete this cancellation might be. 
Moreover, even with a realistically small but finite quark mass, the qq condensate becomes 
of relative size ~ 10% by = 1 GeV^. We therefore conservatively choose Q"^ = 2 GeV^, 
estimating corrections to the lowest order result to be no more than 15%. We thus obtain 
Xy(2.0 GeV^) = (3.52 ± 0.53) ■ lO'^ GeV'^ 



(GeV^) 


«s(g')/vr 




4.0 


+0.101 ±0.013 


-(0.008- 0.025) 


3.0 


+0.109 ±0.015 


-(0.015- 0.044) 


2.0 


+0.126 + 0.020 


-(0.033- 0.099) 


1.0 


+0.175 + 0.037 


-(0.132- 0.396) 



Table 1. Values of corrections to XviQ"^) appearing in Eq. ( 16. 1\) relative to leading order. 
CoefEcients have been suppressed for simplicity. 
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The isospin factor n/ in Eq. (f4.2|) for the pion case is set to unity; in previous cases, more 
than one isospin channel could couple to the current in the same dispersion relation 
so that n/ > 1, thus strengthening the expression (p.9|). In the tt case this is no longer 
possible. As we have pointed out, tt" pairs cannot couple to a vector current due to charge 
conjugation. Moreover the electromagnetic current possesses both J = and 1 amplitudes, 
so an improvement through ri/ with other contributions to the dispersive bound would first 
require disentangling these amplitudes. 

The dispersive bounds may also be strengthened by the inclusion of perturbative QCD 
information along the cut in the deep Minkowski region t — > +cxd, as in However, we 
opt not to do so in this work, for in this approach one must not only select a point t = 
for the onset of this region, but employ the perturbative expression over the whole interval 
(t*,+cx)) for inclusion in the dispersion integral. The uncertainties are therefore those of 
the integrated perturbative result, rather than those of just one point in the Euclidean case. 
For our purposes here, we prefer the more minimal approach of using only deep Euclidean 
QCD calculations. 

We select the parameter to = 0. This natural choice means that 2; = occurs where 
current conservation normalizes the form factor F{t = 0) to unity, the charge of the pion. It 
further means that the form factor near this point is well determined by the first few 6„'s, 
owing to the geometric convergence of J^n^nZ"'- Of course, such convergence is contingent 
upon the small size of |6„|, which is not so obvious as that of the \an\ [compare ( p.lO|) and 
(^.5|)] but is empirically true in our fits. 



Explicitly, the normalization of the form factor and the pion charge radius (r^) 
6 {dF/dq"^)] are given by 



\F{t 



[^1 + Q2/4m2 + 1 



(6.2) 



The experimental data for the form factor in the timelike region is collected from a num- 
ber of sources 0, and each data point is regarded as having significance entirely determined 
by its stated uncertainty. This set contains many dozens of points stretching from = 0.1 
to almost 10 GeV^ and covers decades of experiment, although the most recent published 
measurements are already over ten years old.0 To be precise, our sample consists of 145 
timelike data points, which gives the absolute upper limit of the number of parameters 
bji one may hope to extract from the data; more on this in a moment. 



Recent spacelike results from E93-021 at Jefferson Lab are not yet in print and thus have not 
included in the analysis. 
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B. The Fit 



We minimize x with regard to the unknown parameters 6„ subject to the sum rule 



by using the Levenberg-Marquardt method [11]. Naturally, the topology described by a 



function with scores of parameters utterly escapes intuitive expectations about where the 
minimum might lie, so several choices of parameters for the initial iteration are selected to 
test the proper convergence of the algorithm to the global minimum. Typically, on each given 
run the program finds the same Xmin within a few tenths of a percent. Uncertainties on 
the parameters are estimated by means of the covariance matrix, despite the fact that ( ^.4|) 
is quadratic rather than linear in the parameters Indeed, the two sets {±6„} trivially 
give the same fit, but we select 60 > to guarantee the positivity of the form factor at 
t = to = 0. The program is permitted to wander freely in parameter space without 
imposing the constraint ( |5.4| ); the degree of saturation of this bound, which expresses quark- 
hadron duality, is computed at the end of the fit and is observed to satisfy the bound in 
all cases. The 15% uncertainty in the perturbative value of appears as a systematic 
uncertainty in the parameters 6„ since oc so that the dispersive bound ( ^.5|) 

has a linear uncertainty in Ax^ , meaning that the 1 on the r.h.s. is to be replaced by 1.15. 

The figure of merit in these fits is P^r degree of freedom (d.o.f.), where the number 
d.o.f. is the number of data points minus the number of fit parameters N. Obviously the 
average variance of the fit value of the form factor from data decreases as is increased, 
but the d.o.f. denominator eventually compensates for this advantage, and so there exists a 
particular value of such that x^/d.o.f. achieves a minimum. For our sample of 145 data 
points, this occurs for about 60 parameters 6„. 

There is another reason to choose substantially smaller than the number of data points. 
If we think of the fit as being essentially a decomposition into Fourier modes (as evidenced by 
the fact that our basis functions are = e™^ on the unit circle), then extracting a number 
of parameters comparable to the number of data points is analogous to probing modes 
with wavelengths as small as the spacing between these points. Spurious high-frequency 
oscillations then appear in the fit, and the same phenomenon is observed in the present 
case. If one fits not to isolated data points but a model with a continuous prediction for the 
form factor, then one is free to expand to arbitrarily high harmonics; the oscillations are an 
effect of finite experimental resolution. 

In fact, these spurious oscillations tend to be compounded with even more dramatic 
effects when we attempt to extrapolate from the data-rich regions around the p peak down to 
t+: The oscillations assume huge proportions, many times larger than the p peak itself, with 
the few data points far below the peak nestling themselves in the minima of these oscillations 
(Fig. 1). The cause of this particular misfortune has a very simple origin: Analyticity in z 
is precisely not the same as analyticity in t, and the map (p.5|) is not conformal at the point 
t = t+, i.e., z = —1. This is obvious from geometric considerations: whereas in t space 
segments of the real axis below and above t+ are parallel, the corresponding segments of the 
unit disk in z are the real axis and the two halves of the unit circle, which make 90° angles 
in antiparallel directions at z = —1. Indeed, the Jacobian of the map is given by 

dt 4(t+-to)(l + ^)' ^ ' 
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meaning that any parameterization in the variable z has problems as t t+. And yet this 
is a necessary feature of the dispersive analyticity constraints: It is z, and not t, that is the 
natural kinematic variable for describing the analyticity of the form factor. Additional phys- 
ical requirements or a greater density of data are necessary in order to impose smoothness 
on the form factor in the region near t+. 




0.5 1 1.5 2 

t (GeV^) 



FIG. 1. Best fit oi the timelike tt elastic form factor data \F{t)\'^ to the parameterization 
of ( |5.^ j with parameters {bo, bi, . . . , feeo}- The range < t < 2 GeV^ is shown, although 
some data exists out to 10 GeV^, in order to emphasize the p peak. Note especially the 
difhculty of the fit in accommodating the u shoulder, oscillations of the fit below the peak, 
and the tendency of the fit to become unhinged for t near threshold where there are gaps in 
the data. 
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C. Results 



Our best fit to the form factor over the timehke region appears in Fig. 1. As stated 
above, the fit uses = 61 parameters {bo,bi, . . . ,fe6o} and has a xLin/^-^-^- 3.20. This 
large value, despite the visual goodness of the fit in the figure, can be explained by the 
mild inconsistency of data from numerous different experiments, as well as the famous "u; 
shoulder" due to the presence of destructive interference between the narrow u and broad 
p peaks. The shoulder represents a sudden, nearly discontinuous change in the form factor, 
which is difficult to simulate using only the lower harmonics. Indeed, eliminating the 15 
data points in the immediate vicinity of this shoulder from the fit diminishes x^j^/d.o.f. 
by more than 1/3; of course, the oscillations in the fit below the p peak remain. A third 
possible explanation is that is large because ^ 60 might be required to fit the data 
adequately. However, suppose one considers not Xmin/^-o.f. but rather Xmin / dO'tum and 
fit to 200 parameters (approaching the limitations of modern workstations). Ignoring the 
numerous spurious oscillations produced by this stretch, we find that Xmin/ datum decreases 
from 1.85 to 0.75, i.e., much slower than linearly with N. 

The degree of saturation in the best fit, namely, the fit value for the l.h.s. of ( p.4| ) 
divided by the loosened bound of 1.15, turns out to be 0.466, meaning that charged pions 
and resonances coupled to them account for about half the dual perturbative result at 
Q2 = 2GeV^ 

Using Eq. ( |6.2| ), the values for the spacelike parameters extracted from the = 61 fit 
are given by 

= 0) I = 2.56 ± 2.00, (r^) = 68.2 ± 88.3 GeV"^. (6.4) 
On the surface, this is not a very impressive set of results, considering the values extracted 



directly from the spacelike data |jTO| 



|F(t = 0)1 = 0.995 ±0.002, (r^) = 10.26 ± 0.26 GeV-^ (6.5) 

What has gone wrong? To sharpen this complaint, we observe that even after fitting to 
a full 60 parameters, one still cannot do a very good job extrapolating the timelike data 
into the spacelike region. After all, the naive behavior of the data plotted in t seems to 
suggest that the lower edge of the p peak extrapolates smoothly into the spacelike region 
to give much better figures for the normalization and slope of the form factor at t = 0. 
However, as we have pointed out above, the dispersive bound uses analyticity in not 
t. Models based upon the expected theoretical shape of the p peak or chiral perturbation 
theory implicitly assume that no peculiar behavior afflicts the form factor near For the 
dispersive bounds to do the same, one would have to absorb the {l + z)~^ factor in ( |6.3D into 
the parameterization in order to make dF/dt finite at However, since we have insisted 
upon z as the natural variable of analyticity, for the remainder of this work we maintain 
that such a modification takes us too far from the original motivation of rigorous, minimal 



^In fact, one might expect a discontinuity in dF/dt at indicating the threshold of absorptive 
processes entering through ImF. 

17 



bounds based on quark-hadron duality and analyticity. Certainly such modifications are 
straightforward to implement, but we forgo them for now for the sake of minimality of 
assumptions. 

The next issue is how one can believe the extrapolation from the timelike to the spacelike 
region if indeed the best fit is pathological in the neighborhood of Again, the answer is 
that it is z and not t that is the relevant variable of analyticity. An extrapolation along the 
real t axis, expressed in z coordinates, consists of following the contour of the unit circle in 
z until reaching [z = —1), and then moving along the real z axis to t = (2; = 0). Such a 
path takes us straight through the eye of the storm at t+, and should be avoided if possible. 
In fact, since the form factor is analytic on the entire unit disk in 2;, one may choose a 
more direct and less contentious path; for example, start at the point on the circle where 
data is plentiful, such as tp = tim^^ (corresponding to = 42.5°), and proceed directly 
along a radius to t = 0. The fit for |Fp along this contour is presented in Fig. 2. This 
contour expressed in complex t space resembles a cardioid with its outer edge at tp and its 
cusp at i = 0. Since it does not pass near it does not exhibit strong oscillations, as is 
apparent from the figure. Clearly, analyticity in z has no trouble with this region, and so 
the extrapolated results have meaning. 
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So our fit values are meaningful, and the numbers extracted for \F{t = 0)| and (r^) are 
certainly consistent with those from the spacelike fit, but why are the uncertainties so large? 
Again, the interpretation becomes clear in z space. The position of points z{t, to) in the 
unit disk is completely determined by the ratios t/t+ and to/t^, as is clear from (P75|). With 
our choice to = 0, the timelike data becomes compressed almost entirely into the half-circle 
Re z > 0: See Fig. 3. As mentioned above, the p peak data is clustered about 9 ~ 7r/4, 
while the lowest measured points are at t ^ 0.1 GeV^ —>■ 9 27r/3. This is a direct result 
of tp ^ (t+ — to)- If one takes to < 0, the extrapolation to to should be much more precise. 
Unfortunately, to = is exceptionally convenient for the extraction of \F{t = 0)| and (r^), as 
we have discussed above. We expect these problems would become much less severe if more 
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near-threshold pion production data were measured. That this is true may be confirmed by 
performing the fit after including additional Monte Carlo data points in this region. The 
extrapolation uncertainties would also decrease substantially if we considered form factors 
where much more data naturally falls just above threshold; this is precisely the situation for 
the K form factor, where the peak occurs very close to the threshold Am\, although in 
that case one faces the prospect of depleted data for 9 < 7r/2. 




o 



FIG. 3. A density plot of timelike data on the unit circle in z. Note the preponderance of 
data near the p peak (9 — 42.5°j and the paucity of data near threshold t+ (9 — 1S0°), or 
even beyond 9 — 90°. 

The distribution of the parameters bn in magnitude is an interesting test of the conver- 
gence of the series J2n i>nz"'- Recall that we no longer have specific mathematical information 
on the rate of the convergence of this series. In our best fit with iV = 61 we find only 3 
parameters with > 0.05, 12 with 0.02 < < 0.05, 11 with 0.01 < \hn\ < 0.02, and 
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35 with \bn\ < 0.01. The largest single parameter is 60, and the largest 8 parameters all 
occur in {bo, . . . , 610}. Although the series is not monotonic, it appears to exhibit behavior 
consistent with fairly rapid convergence in n. 

It is also possible to perform the fit to data fixing the normalization \F{t = 0)| = 1 by 
choosing 60 via Eq. (|6.2|). In a fit with \F{t = 0)| = 1 using spacelike data, [|10| gives 

(r^) = +11.07 ±0.26 GeV"^ (6.6) 

whereas our fit with timelike data gives 

(r^) = -13.3 ± 27.9 GeV"^. (6.7) 

Again, our fit value is consistent with the other determination. As expected, the uncertainty 
in (r^) decreases once \F{t = 0)| is fixed, since an additional constraint has been placed on 
the set of allowed analytic functions that satisfied the original fit. The value of Xmm/^-'^-^- 
for this fit is again 3.20, and the level of saturation is now 0.52. 

If enough data existed that the uncertainties on the parameters bn were much smaller, one 
could extrapolate throughout the entire spacelike region t G (— cxd,0), which is z G (0, +1). 
The prediction for |Fp for our best fit values turns out to fall smoothly and monotonically to 
zero as t — > cx), and we expect this paradigm to persist in similar fits with more constraints. 
Because of the large uncertainties in our extrapolations, we decline to exhibit just the central 
value fit; however, we expect this paradigm to persist in similar fits with more restrictive 
constraints. 

Finally, the model-independent parameterization of (|5.3|) may be used directly in a fit to 
spacelike data. This application is entirely analogous to the semileptonic fits of 0-0], except 
that we pick to = and do not here seek to minimize the "truncation error" of between the 
fit and true form factor at points with t <^0: After all, we are interested in moments about 
t = 0. A collection of 45 data points in the spacelike region taken from [^] may be fit with 
just {bo, bi, 62} to give Xmin/d-o.f. of only 0.97 and a saturation of the dispersion relation of 
only 0.018 (as expected, fewer parameters lead to less saturation of the dispersive bound). 
From the parameter fit values we obtain 

|F(t = 0)1 = 0.997 ±0.005, (r^) = 11.28 ± 1.86 GeV-^ (6.8) 

or taking \F{t = 0)| = 1, 

(r^) = 12.33 ± 0.51 GeV-^ (6.9) 

with x^jj^/d.o.f. = 0.98 and a degree of saturation 0.022. As a parting shot we point out 
that this model-independent determination of (r^) differs from that obtained by a naive 
linear extrapolation in t ( |6.6| ) by over 2a. The difference lies in the fact that while the usual 
determination assumes F linear in t near t = 0, the model-independent form factor has a 
natural curvature associated with the factor l/(j){z), which in turn arose from analyticity 
constraints. The larger size of the uncertainty in our determination takes into account 
theoretical uncertainties in the shape of the form factor, which of course are absent when 
one assumes linearity in t. 
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VII. CONCLUSIONS 



We have explored some of the constraints that timehke data and dispersive constraints 
place on the shape of form factors, particularly in the spacelike region. Much can be learned 
from little more than the analytic structure of the two-point function to which the form 
factor contributes and a few mild physical assumptions from finiteness. The sum rules and 
improved parameterization derived above arise directly from these minimal considerations. 

Nevertheless, the extrapolation of a finite amount of timelike data to the spacelike region 
is seen to lead to huge uncertainties in the case of the pion elastic form factor F{t). This 
means that the bare assumptions of analyticity and the applicability of QCD in the deep 
Euclidean region are not enough to obtain what seems to be a natural extrapolation of 
the p peak to the spacelike region. We have discussed the dangers of extrapolating in the 
kinematic variable t rather than the natural kinematic variable of analyticity, z: Namely, 
one must pass around a branch point t+ — 4m^ or z — —1, near which the smooth behavior 
of the form factor may be compromised. In terms of z, the assumption of a resonant pole in 
the timclike region places extremely tight constraints on the allowed parameter space of 6„. 

On the other hand, the z contour from timelike to spacelike data may be distorted to 
avoid the troublesome point z = —1. The large uncertainties in our fit are seen to arise from 
the density and placement of the data set in z; basically, the problem is that most of the 
data, clustered near the p peak, is far above the pair production threshold t+. Since this 
is not the case for the K form factor, the extrapolation may be much more reliable in that 
case. 

This is not to indicate that the technique is hopelessly weak or intrinsically flawed, 
however. Rather, we have endeavored to exhibit what is obtained from the absolute minimum 
of assumptions. We have also suggested at various points possible improvements to the 
program presented above. It is useful to collect them here for the convenience of the reader. 
First, the calculation of the perturbative quantity X^(Q^) be improved by a careful 
analysis of multiloop effects and condensates. We have evaluated the perturbative result 
only at the single value = 2 GeV^. Certainly the behavior of the dispersion relation as 
a function of Q^, which is equivalent to a moment analysis, provides additional constraints. 

Second, the stringency of the dispersive bound is determined by both x^iQ'^) fhc level 
of saturation by all contributing states excepting those being probed by the parameterization. 
In the present case, this means any states like KK that receive contributions from the 
electromagnetic current but have not already been counted in n^7i~ production (such as 
the p) . A calculation or estimation of this contribution effectively decreases the dispersive 
saturation limit allotted to the pion form factor. Also falling into this category is the 
deep Minkowski region t ^ of the pion form factor, for which data has not yet been 
collected. Here, one may either use its integrated weight to increase the stringency of the 
dispersive bound, or retain its functional form to obtain another sum rule on the form factor 
parameterization . 

The latter suggestion leads us to the third possible improvement, that conditions on 
the shape of the form factor lead to constraints on the form of the parameterization. We 
have endeavored to make as few assumptions as possible along these lines, lest the model- 
independency of the parameterization is lost. Nevertheless, modelers are free to start with 
the dispersive bounds as a starting point for building form factors that automatically satisfy 
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nontrivial QCD constraints. 

Another possibility along these lines not yet mentioned is the choice of the series used 
to parameterize the unknown physics of the form factor. We have always used a simple 
Taylor series in to express the analyticity of the form factor since the set {z''} is a 
complete and orthonormal basis on the unit circle for analytic functions. Although this basis 
proved immensely useful (indeed, optimal) for spacelike fits, perhaps a better basis exists for 
accommodating the relatively rapid variations in the form factor due to structure in resonant 
peaks, such as the lu shoulder. However, the constraints on this new parameterization are 
manifold: If it is to be adequate for both timehke and spacehke fits, it must admit a Taylor 
series rapidly convergent in as 2; — > 0, all of its Taylor coefficients must be real (so that 
the form factor is real in the spacelike region, where z is real), but it cannot converge too 
quickly in 2;" as |2;| ^ 1, or else our polynomial fits in z"' would already have been adequate. 
Moreover, since the proposed functions are analytic in z and not z on the unit disk, natural 
choices such as cos n9 and sin n9 are forbidden. 

Clearly, many aspects of this program yet remain to be fully explored, but the variety 
of directions for improvements are quite astounding. Dispersive techniques possess the rare 
ability to encapsulate a great deal of nontrivial physics with true simplicity and elegance. 
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